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simulation  and  experimental  investigation  is  used  to  achieve 
a basic  understanding  o£  diffusion  and  optical  wave 
propagation  in  glass.  Normalized  propagation  constants  and 


asymmetrical  graded  index  planar  waveguides  and  graded  i 
channel  waveguides  with  arbitrary  index  profiles.  One 


applied  electric  field,  are  solved  numerically  to  determine 
the  nature  of  the  index  profiles  of  the  practical  glass 
waveguides.  The  results  of  diffusion  simulation  are 


profile 


by  a new  and  simple  technique — scanning  electron  microscopy. 

A novel  electrolytic  process  is  developed  for  releasing 
silver  ions  in  molten  sodium  nitrate  and  a Galvanic  cell 


silver  ion  concentration. 


parameters  are  derived  and  used  subsequently  to  formulate  a 
systematic  procedure  for  fabricating  single  mode  waveguides 
entirely  in  terms  of  the  process  parameters.  Finally, 
symmetric,  buried  waveguides  with  a fiber-like  index  profile 
are  fabricated  using  a two-step  diffusion  process. 


Optical  Communication  Systems 


edium  for  a number  of  present  and  future 
offer  many  significant  advantages  compared  to  other 


waveguides).  T 


and  circular  and  rectangular  metallic 
include  cost  effective  long  distance 


high  bandwidth  capabilities,  small 
degree  of  intercept  security  and 


size  and  weight,  high 
dielectric  isolation 


The  foregoing  desirable  features  of  optical  fibers  have 
strongly  stimulated  efforts  in  associated  technologies  such 
as  fiber  cabling,  long  life  solid  state  sources,  high 
performance  detectors,  integrated  optical  couplers  and 
modulators.  These  technological  advances  have  propelled  the 

1970s  into  a major  industry  today.  The  ongoing  commercial 


installations.  field  trials  and  ] 
n future  telecommunication  a 


research 
l computer  systems. 


The  basic  functions  of  an  optical  communication  system 

figure  1.1.  The  system  is  used  for  transmitting  analog 
signals,  e.g..  speech  or  video  signals  or  digital  signals, 

coded  using  some  modulation  scheme  such 
Modulation  (PCM)  or  Delta  Modulation  (dm). 

pulses.  The  optical  fiber  is  used 
pulses  to  the  receiving  end  where  a photodiode  is  used  to 


3 modulator 


dispersion,  which 


can  successively  degrade  the 
fiber.  The  impurities  in  the 


resulting  in  the  possibility  of  causing  inter-symbol 
interference.  This  necessitates  the  use  of  repeaters  at 

regular  intervals  along  the  fiber,  where  the  signal  is 


boxes  A and  B shown  in  figure  1.1  represent 
a class  of  devices  called  the  integrated  optical  devices 
which  have  been  recently  introduced  in  the  conventional 

signal  processing  components,  in  the  sense  their  inputs  and 
outputs  are  optical  signals.  They  can  perform  a variety  of 

electrooptic  signal  processing  applications,  e.g.,  switches 
and  modulators,  an  electrical  signal  containing  the 

output.  Introduction  of  these  devices  has  considerably 
improved  the  performance  of  optical  communication  systems, 
as  illustrated  in  figure  1.2.  with  the  possibility  of 
multiplexing  and  demultiplexing  several  wavelengths  on  the 


carrier  optical  communication  system  can  be  multiplied  by 

Wavelength  Division  Multiplexing  (WDM) . Another  area  where 

fibers,  also  shown  in  figure  1.2.  By  using  two  different 
wavelengths  in  opposite  directions,  a single  fiber  can  carry 


fiber-optics  technology 


cabling 


Excellent  review  papers  covering  all  these  aspects  are 
available  in  literature,  for  example,  Hiller  et  al.  [Mi731 , 


1 .2  Optical  Computing 

Researchers  have  been  trying  to  design  an  all-optical 

computer  technology  is  approaching  the  speed  limits  imposed 
by  the  si  or  GaAs  microelectronics  technology,  an  optical 
computer  is  expected  to  improve  the  performance  by  orders  of 


system  cycle  time  and  pulse  widths  shrink,  the  bandwidth 
needed  to  preserve  the  rising  and  falling  edges  of  the 
pulsed  signals  increases.  This  makes  it  necessary  to  have 


possible.  For  example, 
length  of  6 inches. 


Optical  components,  in  principle,  provide  a possible 
solution  to  this  problem  because  of  their  large  bandwidth 
and  non-interfering  propagation  [Hu84] . Difficulties  in 
realizing  all  optical  logic  devices,  which  are  the  building 


problem. 


Recently,  however,  a feasible  design  for  optical  AND  and  XOR 


picoseconds 


gates  with  propagation  delays  as  low  as  10 

optical  bistability  have  also  been  proposed  [Sm81,  Ta81]. 

optical  non-linearities  of  dielectric  (electrooptic) 

such  as  waveguides  and  directional  couplers  are  likely  to 
become  the  basic  components  for  the  hardware  of  the  future 
optical  computers.  Architectures  for  such  high  speed 
optical  computers  have  already  been  proposed  [Ta79b] . Their 

is  capable  of  performing  addition  and  multiplication  without 
the  necessity  for  carry  operation.  Calculations  can 
therefore  be  performed  in  a fully  parallel  manner  utilizing 

Besides  telecommunications,  integrated  optical  devices  are 
expected  to  find  a large  number  of  applications  in  this  new 
and  exciting  area. 


progress  over  the 
primarily  on  the  £ 

configurations,  < 


of  the  Progress  in  Integrated  Optics 
i in  Integrated  Optics  has  continued  to 


optical  materials. 


sm  into  appropriate 
ridge  and  buried 


channel,  a variety  of 
of  performing  a wide  r< 

miniaturized  optical 
devices 


nge  of  operations  on  light  waves  have 

in  lead  to  extremely  efficient 
circuits  with  low  drive  power 
e light  energy  propagating  in  these 


phenomena  produced  by  high  intensity  fields  such  as 
non-linear  optical  interactions  can  be  observed  at  modest 
input  power  levels.  This  feature  has  added  another  new 

Optics.  As  the  research  is  inter-disciplinary,  involving  a 
number  of  scientific  and  technological  disciplines,  its 


investigators  having  different  backgrounds.  Thus  engineers 
and  scientists  with  interest  in  optics  and  microwaves. 


semiconductors. 


crystal  and  surface  phenomena  have  often  worked  together  to 
explore  new  ways  to  design  and  fabricate  these  devices. 


the  principles  of  dielectric  waveguides — for  example, 
methods  to  couple  energy  into  film  waveguides  (Ti71,  Ko75] , 
and  guided  wave  devices  for  communication  applications 
IA1S11.  The  research  can  be  broadly  divided  into  two 
catagories  {with  the  requisite  overlaps). 


Development 


materials  and  fabrication 


A variety  of  substrate  materials  such  as  polymers. 
Niobate,  Lithium  Tantalate.  Zinc  Oxide.  Ill  - V 


Semiconductors,  etc.  have  been  used  over  the  last  decade. 
The  associated  fabrication  processes  include  ion  implan- 
tation. ion  diffusion,  thermal  outdiffusion,  metal 
indiffusion,  sputter  coating  and  surface  layer  growth. 

Some  of  the  early  films  capable  of  guiding  light  were 
reported  by  Goell  and  Standley.  In  their  process.  Corning 
7059  glass  was  RF  sputtered  on  a microscope  slide  to  produce 
a waveguiding  film  [Go69al . Many  varnishes  and  epoxies  have 
also  been  deposited  on  glass  substrate  to  form  waveguides. 


and  photoresist  made  by  deposition  on  glass  slides  [U1721 . 
Ramaswamy  and  Weber  have  reported  polymer  films  with 
adjustable  refractive  index  [Ra73] . proton  implantation  has 
been  used  by  Schineller  et  al.  to  create  waveguiding  layers 
in  fused  silica  [Sc68) . Kaminow  and  Carruthers  have  used 
outdiffusion  of  Li  atoms  to  form  waveguiding  layers  in 
LiNb03  and  LiTa03  [Ka731.  Indiffusion  of  various  transition 
metal  elements  has  been  used  by  Schmidt  and  Kaminow  to  form 
low  loss  planar  waveguides  in  LiNb03  [Sc741 . 


waveguiding  layers 


silicate  glass — Giallorenzi  et  al.  [ Gi 7 3 J , and  borosilicate 
glass — Izawa  and  Nakagome  [Iz721. 

and  ease  of  fabrication  have  been  the  major  criteria  for 
accepting  or  rejecting  a substrate  material.  From  these 

exchange  in  glass  have  become  the  most  popular  techniques 


The  devices  have  evolved  from  simple  planar  waveguides 
to  complicated  analog  to  digital  converters.  successful 

[Ma73b|,  Standley  and  Ramaswamy  tst74 1 ) was  a major  step  in 
the  development  of  Integrated  Optics  because  it  paved  a way 


stepped  as  reversal  (Kogelnik  and  Schmidt 
interferometers  {Martin  IMa75] ) and 

devices  developed  in  the  late  seventies. 

The  present  thrust  of  research  in 
super  high  speed  electrooptic  modulators  and  fast 

(Taylor  [Ta79a) , Leonberger  et  al.  [Le  82!).  Another  area 
of  active  research  is  device  integration.  The  ability  to 
efficiently  couple  light  between  a single  mode  fiber  and  a 


[Ko76] ) , 
balanced  bridge 


incorporated  into  practical  systems.  Increasing  efforts  in 
reducing  the  coupling  loss  have  resulted  in  the  achievement 

1 = 1.32  urn  has  been  demonstrated  (Ramaswamy  et  al.  IRa82]). 

primary  phase  of  development  in  this  exciting  new  field 
seems  to  be  successfully  completed.  With  the  design  of 

coupling  efficiency  between  these  components  and  optical 
fibers,  the  next  goal  seems  to  be  the  integration  (in  a 

multiplexers,  couplers  and  demultiplexers  on  the  same 
substrate  to  form  functional  sub-systems. 

Figure  1.3  illustrates  a transmitting  end  optical  sub- 
system which  typically  represents  box  A of  figure  1.1.  At 
the  transmitting  end,  an  external  modulator  is  required  if  a 
laser  is  used  as  the  light  source.  The 

ch  exhibit  electrooptic  effect,  e.g.,  LiNbo3.  To  make 
modulator  efficient,  l.e.,  to  have  a small  drive 
kept  as  small  as  possible. 

. strip  waveguides  about  3- 
■ split  the  modulator  output  to 
star  coupler  is  used  as  shown 


Typically  i 


several  parallel  channels. 


Figure  1.4  illustrates 

receving  end.  there  is  no  ne 
passive  substrate  material  s 

compatible  with  a single  mode 

this  sub-system  include 


typical  receiving  end  s 

liber,  thus  eliminating  a n 
Typical  devices  used 
wavelength  filter  and  a s 


coupler  a 


1.4  Project  Definition 

The  majority  of  the  devices  described  in  the  last 
section  have  been  fabricated  on  a LiN'b03  substrate  using 
Titanium  indiffusion.  Since  LiNb03  is  an  electrooptic 
material  the  refractive  index  of  the  guide  and/or  the 
surrounding  region  can  be  altered  by  the  application  of  an 
electric  field.  This  electrooptic  tuning  has  been  extremely 
useful  in  meeting  the  desired  device  specifications.  For 
example,  the  center  wavelength  of  a demultiplexing  filter 

desired  [A178] . Being  a crystalline  material,  the  cost  of 
LiNbOj  is  prohibitive  as  the  growth  of  its  single  crystal 
requires  a carefully  controlled  process.  In  addition,  it 


refractive 


optical  fiber  (~  1.46),  resulting  in  high 
losses  due  to  Fresnel  reflection  at  the 

essential  where  modulation  or  any  other  active  operation, 

of  an  optical  communication  system,  as  there  is  no  need  to 

s glass  could  be  considered.  Refractive 
.me-silicate  glass  is  around  1.5  making  it 
optical  fibers.  In  the  absence  of 
the  "passive"  devices  made  on 
ed  to  compensate  for  the  fabrication 
:t  fabrication  tolerances  and  precise 
needed  to  meet  the  desired  device 


compatible 
electrooptic  effect, 
glass  cannot  be  tun 


specifications. 

In  1972,  Izawa  and  Nakagome  first  established  the 

for  making  graded  index  optical  waveguides  (Iz72].  Since 
then  several  sporadic  attempts  [St77,  Ch80,  Li82]  to  revive 
this  technique  have  occurred;  none  of  which,  however,  has 
resulted  in  reproducible  fabrication  of  single  mode  wave- 
guides. The  process  basically  involves  a bath  containing  a 
mixture  of  molten  AgNOj  and  NaNO^  in  which  soda-lime- 
silicate  glass  slides,  about  1 mm  thick,  are  suspended  for 
Na+-Ag+  exchange.  The  exchange  of  sodium  ions  in  the  glass 


higher  refractive  index  capable  of  guiding  light.  The 
existing  process  lacks  control  on  some  of  its  parameters  and 
no  attempts  are  made  to  improve  the  process  yield. 

is  investigated  with  emphasis  on  the  "reproducibility"  and 

theoretical  part  of  the  research  has  been  directed  first,  to 
the  basic  understanding  of  the  Ag+  diffusion  (exchange) 
process.  This  is  followed  by  the  study  of  electromagnetic 
wave  propagation  in  graded  index  planar  glass  films  so  that 
the  expertise  gained  would,  facilitate  the  design  of  future 
optical  components.  The  experimental  work  is  aimed  at  the 
development  of  a controlled  process  for  reproducible 

the  process  parameters,  such 

(waveguide  width  and  index  profile)  wh 
propagation  characteristics  (number  o 
their  propagation  constants).  Thus,  ii 

necessary  to  find  the  links  relating  the  process,  device  and 

for  the  case  of  planar  guides  formed  on  the  surface  of  glass 
substrates. 


diffusion  process, 
vice  parameters 
s supported  and 


overall  problem  is  divided  ii 


following 


Propagation  analysis 
Diffusion  simulation 


Waveguide  characterization. 

, the  work  carried 
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CHAPTER  2 

PROPAGATION  ANALYSIS 


Propagation  analysis  of  dielectric  waveguide  mainly 
consists  of  determining  the  number  of  modes  supported  by  the 
waveguide,  the  mode  propagation  constants,  and  the  trans- 
verse field  distribution  associated  with  these  modes.  This 

of  integrated  optical  components.  In  end-fire  couplers,  the 
end  faces  of  waveguides  and/or  fibers  are  polished  and 
r without  an  index  matching  fluid  in  between. 

he  two  waveguides.  Thus,  to  achieve  good 
field  profile  in  these  guides  should  have 
d dimensions.  Hence,  to  reduce  the  coupling 

guides  is  essential.  In  directional  coupler-type  devices, 
the  waveguides  are  brought  adjacent  to  each  other  sideways 
and  the  coupling  takes  place  through  the  interaction  of  the 
evanescent  fields  due  to  the  traveling  waves  in  the 
guides.  This  is  also  referred  to  as  “optical  tunneling." 
Complete  energy  transfer  can  occur  only  when  the  propagation 
constants  of  the  modes  are  the  same.  Thus,  determination  of 

device  parameters  are  necessary  for  efficient  device  design. 


The  coupling 

efficiency, 
similar  shape 


19 


Graded  index  structures 
( waveguides  are  formed  in 
suits  in  waveguides  wit 


graded  refractive  index 

1 conditions  manifest  themselves  in  different  index 
profiles.  For  example,  one-dimensional  diffusion  with  an 
Infinite  source  of  dopant  results  in  a complementary  error 


source,  i.e.,  the  dopant  s 
during  diffusion,  yields  a 
For  the  waveguides  formed  bi 
is  essentially  like  a step  j 

practice,  are  considered 
sections. 


i exhaustible  input 
depleted  entirely 
ex  profile  [Bo63] . 
sputter  deposition,  the  profile 

in  detail  in  the  subsequent 


1.  Planar  Step  Index  Waveguides 

2.  Planar  Graded  Index  Asymmetrical  Waveguides 


of  three  layers:  substrate,  film  and  cover  (or  superstrate) . 

used  to  coat  the  waveguide  as  is  the  case  with  polymer  clad 
waveguides.  The  guided  mode  propagation  can  be  understood 

travels  in  a zig-zag  fashion  inside  the  film  undergoing 


Ray  Analysis 


As  the  ray  travels  through  path  ABC  in  figure  2.1,  it 
suffers  phase  shifts  at  both  interfaces  due  to  total  inter- 
nal reflection  as  well  as  a phase  change  while  propagating 


an  integral  multiple  of  2v,  then  the  phase  front  after  total 
internal  reflection  at  A as  it  advances  to  C comes  in  phase 
with  the  phase  front  which  is  reflected  at  C,  resulting  in  a 
constructive  interference.  Under  these  conditions,  the  film 
supports  a ray  corresponding  to  a specific  angle  0.  For  a 
given  set  of  film  parameters,  the  rays  propagate  at  specific 


condition  can  be 


where  m is  an  integer  and  the  phase  shifts  es  and  eQ  can  be 
expressed  in  terms  of  film  parameters  ITi71,Ko791.  For  TE 
modes  with  the  electric  field  polarization  parallel  to  the 


arguments 


2 and  <V"c ' 2 


"£'  n„,  d,  and  XQ  v 


n£2  - nb2 


■-3*^ 


• nb  and  "c  <(  ne 


/(l-b)^  = tan_1(b/U-b)|^  + tan-1  ( <b+a)/(  1-b)  + mt 


which  are  similar  to  8-w  diagrams  of 

that  they  represent  normalized  propagation 
TM  modes.  This  makes  them  helpful  in 
waveguide  design. 


the  rectangular 


2.2  Planar  Asymmetrical  Graded  Index  Waveguide 
The  model  for  a planar,  asymmetrical  graded  index 
waveguide  is  shown  in  figure  2.3.  The  index  profile  of  a 
graded  index  waveguide  occupying  the  half  space  {x  > 0)  is 


(2.9) 


nb  = bulk  refractive  index 
ns  9 surface  refractive  index 

The  index  profiles  for  various  possible  functions  f(x)  are 
plotted  in  figure  2.4. 

profile f (x) 

Complimentary  error  erfc(x) 


exponential 


depends 


fabrication  process. 


2.2.1  Normalization-Dispersion  Equation 

done  almost  on  the  same  lines  as  the  step  index  except  that 


internal  reflections  where  as  in  graded  index  guides  the 
rays  follow  curved  paths  with  a total  internal  reflection 
occurring  at  the  cover-film  interface  and  a turning  point 
occurring  inside  the  graded  index  film  as  shown  in 


propagation  constant  - 


case  is  given  by  the  integral 


(2.10) 


This  is  also  referred  to  as  the  dispersion  (or  characteris- 

the  turning  point  xt.  Hocker  and  Burns  [Ho75]  have  shown 
that  = -x/2  by  considering  the  phase  change  in  the  limit 
0(x)  ♦ x/2.  The  phase  shift  *c,  upon  total  internal 

reflection  at  the  cover-guide  interface,  is  given  by 
equation  (2.3). 

The  dispersion  equation  (2.10)  can  be  normalized  as 


The  following  parameters  are  defined: 
x'  = x/d  = normalized  depth 
n(x')  = nb  + in  f(x') 


The  same  definitions  for 


an  be  used  {equation  (2.4) 
changed  to  ns,  represent- 


V * K0d  (2nb 


Vi 


For  a highly  asymmetrical  case.  nc  <<  ns.  Also  for 

The  rays  in  these  waveguides  are  almost  on  axis,  so 
that  0 - ir/2.  Thus,  sin  6-1  and  cos  6-0.  Substituting 
these  expressions  in  equation  (2.3),  we  get  4 * 

integral  equation  (2.10)  can  be  converted  to 

V / b |f(x'|  - b)  Vidx1  = (m  + 3/4)  » (2.11) 


At  the  turning  point  xt,  th 
direction  of  propagation  (figure  i 


n(x')  - N . 

n(x-)  - n,  + 4n  f(x') 

n2(x')  = *2  = n2  * Un  £<*£)  nfa  , 


(2.14) 


transverse  Helmholtz  equation  [T1711  , 


Depending  upon  the  index  profile  nix),  this  second- 

different  techniques.  Con we 11  has  considered  an  exponential 
profile  and  has  obtained  a closed  form  analytical  solution 
in  terms  of  Bessel  functions  ICo74] . Marcuse  has  approxi- 

segments  and  has  shown  that  the  mode  profile  can  be 

functions  [Ma73a] . For  a quadratic  index  profile,  this 
equation  takes  the  form  of  one-dimensional  Schroedinger 
equation  for  a quadratic  potential  well  and  the  solution  is 
similar  to  that  of  a quantum  mechanical  harmonic  oscillator 
I Ko79]  . The  closed  form  solution  is  expressed  in  the  form 
of  Hermite-Gaussian  functions.  Smithgall  and  Dabby  have 
used  numerical  techniques  to  solve  the  differential  equation 

The  analytical  techniques  provide  closed-form  solutions 

nential  but  cannot  be  easily  extended  to  more  useful 
Gaussian  and  erfc  profiles.  Some  of  the  closed-form 
solution,  e.g.,  for  exponential  profile,  involve  Bessel 


functions  of  arbitrary  order,  Jj>7(x),  J3 ^(x),  etc.,  which 
are  difficult  to  evalute  and  there  are  no  tables  readily 
available.  The  numerical  techniques,  on  the  other  hand,  are 
fairly  straightfoward,  but  the  solution  applies  to  specific 
cases  and  recomputation  is  necessary  each  time  the  waveguide 


differential  equation  (equation 


then  the  solution  bee 
it  is  applicable  to 
index  profile  [Ra83) . 


derivatives 


3 application  of  the  numerical  method, 
>mes  fairly  general  in  the  sense  that 
all  waveguides  with  that  particular 
By  changing  the  variable  from  x to  x1 
usual  definitions  for  the  normalized 


expressed 


% £(*) 


4n)^ 


+ V2[f(x')  - b]  f 


boundary  condition 
condition,  a unit 


eld  can  be  assumed  at  the  cover-guide 
is  of  generality,  because  the  solution 
linear  differential  equation  is  unique 
multiplier.  (This  means  if  Ey(x')  is  a 

Thus,  the  second  boundary  condition  is 


E(  0)  = 


value  was  solved  using  a fourth  order  Runge-Kutta  method. 
Since  for  quadratic  index  profile,  the  solution  exists  in 
closed  form  as  Kermite-Gaussian  functions  (Ko79J,  this 
profile  was  used  to  verify  the  results  of  our  numerical 
method.  The  numerical  results  were  in  excellent  agreement 

of  the  method.  It  was  then  extended  to  the  more  useful 

Figure  2.6  shows  the  b-v  characteristics  for  the  first  three 
modes  of  the  Gaussian  profile.  Figure  2.7  shows  the 

higher  V number,  better  mode  confinement  i 
peak  of  the  modal  amplitude  profile  shifts  toward  the  cov* 

profiles  for  m = 0,  1,  and  2 for  a multimoded  waveguide  w: 
a large  V number,  V - IB.  The  modes  in  the  plot  ; 


s expected. 


normalized  to  ci 


and  thus,  most  of  the  energy  is  confined  to  the  region 
further  away  from  the  surface.  This  confirms  the 

to  the  higher  order  modes,  shows  higher  propagation  loss  in 

prone  to  surface  scattering  loss  [Ra75] . A similar  e 
was  carried  out  for  the  complementary  error  function  i 
profile  and  the  results  are  plotted  in  figures  2.9,  2. 1C 


Since  the  mode  profiles  are 
to  define  a normalized  modal  wid 
distance  between  two  points  of 
intensity  falls  to  1/e  of  its 
effective  diffusion  depth.  Si 
number,  it  is  plotted  in  figures 


normalized,  i 


I 2.13  as  a function 

a Gaussian  index  profile  waveguide  is  single 
used  to  select  t 


diffusion  depth 


regime 


efficient 


end-fire  coupling  between  the  fiber  and  the  waveguide. 

To  understand  how  index  profile  affects  the  field 
profile,  several  waveguides  with  different  index  profiles 


equal  to  8.  The  field  profiles  of  tl 
calculated  numerically  and  plotted  a 
It  is  seen  that  though  the  step  i 


markedly  different 
guides  essentially  h 
fore,  essentially  ha 


profile,  all  other  graded  index 

the  same  normalized  effective  modal 
be  understood  by  considering  the 

number,  the  solution  to  this  equation  depends  on  the 
coefficient  (f(x')  - b).  Figures  2.4  and  2.5  indicate  that 
for  profiles  with  higher  values  of  refractive  index  f{x')  at 
any  x’  also  have  a larger  b for  given  V so  that  the  differ- 
ence (f(x')  - b)  is  approximately  constant  for  all  the 

profiles.  This  result  has  significant  implications  in  the 


of  the  index  profile,  therefore,  does  not  seem  to  be  as 
important  in  coupler  design  as  is  maintaining  the  required  V 
number.  As  pointed  out  earlier,  the  profiles  are  determined 


by  the  diffusion  process.  In  practice,  it  is  very  difficult 
to  maintain  either  an  infinite  source  of  dopants  or  to 

diffusion.  Thus,  the  profile  may  be  somewhere  between 

Gaussian  and  erfc.  The  above  result  indicates  that 
achievement  of  a specific  index  profile  is  not  absolutely 
necessary  for  single  mode  applications.  However,  the 
selection  of  V number  in  conjunction  with  the  diffusion 

the  diffusion  process,  is  extremely  important. 


2.3  Planar  Graded  Index  Buried  Waveguides 
Buried  waveguides  fabricated  in  glass  by  a two-step 
diffusion  process  have  a skewed  index  profile  (Ch80]  which 
can  be  approximated  by  two  linear  segments  as  shown  in 

For  propagation  analysis,  the  following  parameters  are 


' ■ x/d  - normalized  d 


k refractive  i 


waveguide  depth 


refractive  index 


The  usual  definitions  for  V and  b parameters  are 
considered  (equation  2.4  and  2.5),  except  that  the  film 
refractive  index  ng  is  replaced  by  the  peak  refractive  index 

As  seen  from  the  ray  picture  in  figure  2.15,  there  are 
two  turning  points  for  a ray  corresponding  to  a well  guided 

written  as  a modification  of  equation  (2.11)  as 


phase  shifts 


turning  point  x^f  x£2  (figure  2.15), 
as  described  in  the  previous  section. 
This  equation  can  be  split  as 


:he  normalized 


5 linear  functions. 


f,(x')  - (1  - x’)/U  - r 


e turning  points,  as  shown  in  the  equation  (2.1 


equation 


waveguides 

circular 


reported  in  the  literature.  Gc 
numerical  analysis  tor  a 


step  index 


channel  immersed  in 
[Go69b] . Marcatili 

ignoring  the  fields  i 

Hocker  and  Burns  (Ho77] 
pichot's  method  based  ( 
and  Matsuhara  [Ma73c] 
method. 


. channel 


infinite  medium  of  ] 
analyzed  the  step  in< 
refractive  indej 
i corner  region  [Ma69] . 
includes  Effective  Inde 
feh's  Finite  Element  method  IYe751 , 
vector  Integral  equations  [Pi82], 
and  Taylor's  [Ta76I  Variational 


x method  c 


All  these  methods  suffer  from  the  lack  of  generality  in 

specific  waveguides  and  extensive  recomputation  is  necessary 
if  the  waveguide  parameters  are  changed.  Normalization  of 
the  scalar  Helmholtz  equation  prior  to  its  numerical 
solution  is  expected  to  make  the  solutions  fairly  general. 

A new  method  called  a "Variational-Finite  difference" 

waveguides  with  arbitrary  index  profiles  and  aspect  ratios 
fLa84a] . The  method  has  the  following  advantages. 

1.  It  can  accept  any  index  profile,  even  if  it  is 
available  at  discrete  points.  The  method, 
therefore,  can  directly  use  the  results  of 


solution  of  the  two-dimensional  diffusion  equation. 
It  accepts  various  waveguide  aspect  ratios  without 
any  additional  reformulation.  This  feature  is 
important  because  the  aspect  ratio  of  a practical 
waveguide  depends  on  both  the  diffusion  depth  and 

lateral  diffusion.  Thus,  the  aspect  ratio  usually 


3.  The  method  solves  the  normalized  Helmholtz 
waveguides  as  long  as  the  index  profile  remains  the 


The  dielectric  channel  waveguide  under  consideration  is 

the  width  w and  the  depth  d should  be  replaced  by  the 
effective  width  and  the  effective  depth  respectively. 

The  index  profile  can  be  expressed  as 


where  f(x,y),  as  before*  is  the  normalized  index  profile 
function  and  has  a value  between  0 and  1. 

For  the  substrate  region,  n(x,y)  ■ nb,  so  f(x,y)  ■ 0. 

refractive  index. 


For  small  index  difference  (An  <<  nb) , the  modes 
supported  by  this  type  of  waveguide  can  be  divided  into  two 
groups: 

1.  Quasi-TE  modes  in  which  the  main  field  components 


components 


expressed  as 


where  the  field  amplitude  E<x,y)  satisfies  the  two- 


The  following  normalized  space  variables  are  defined: 


w = normalized  w 


y/d  - normalized  depth 


• ' 'J  «*  li*.  &■ 


(x-,y: 


F - V2£<x\yM  E - (^|)2  - 


♦ - / / <V2£(x\y')E2  - -i  (- 


functionals 


/„/  [V2f(x',y')E2  - -5  (jIt)  - (jjt) 
/ / B2dx'dy' 


♦ = I - <v2b>  J 


Comparing  equations  (2.3$)  and  (2.41). 


stationary  value  of  I/J. 

To  proceed  furtherff  it  is  necessary 


finding  the 


distribution  E(x'.y')  and  then  use  the  variational  approach 
to  find  the  undetermined  coefficients  in  the  functional 


Matsuhara  [Ma73c]  and  Taylor  (Ta76)  have  expressed 
E(x.y)  as  a linear  combination  of  a finite  number  of 
parabolic  cylinder  functions 

E(x,y)  =1  a1.l2»?nilj)  D1(x/{)Dj ( y-yQ)/n]  (2.43) 


nd  j respectively. 

ake  the  numerical  realization  easier,  a new  form  of 
is  chosen  using  the  finite-difference  approach.  As 

attains  a negligibly  small  value 

assumption  is  justified.  Obviously  for  modes  closer  to 
iff,  the  distance  needs  to  be  increased  much  further. 


1 guided  modes  t 


E(x' ,y' ) is  defined  at  t 
grid  point  values  (F06O] . 
square  is  written  as 


vn  inside  the  box  and  t 
i grid  points.  Inside  t 

In  figure  2.19,  a typical  square 
the  field  distribution  for  this 


Since  this  form  of  interpolation  is  linear  on  each  link 
of  the  net,  it  provides  E{x',y')  which  is  continuous  on  the 
whole  region.  If  we  select  E(x',yM  to  vanish  at  nodes  on 
the  boundary  PORS,  then  it  automatically  satisfies  the 
boundary  condition. 


The  integral  functions  X and  J defined  by  equations 


III  E2(x',y' jdx'dy'  (2.45) 


is  a single  square  w 
3 single  square. 


Thus,  a vector  U is  defined  whose  entries  are  the 
values  of  E(x',y')  at  the  grid  points;  then  the  summation  in 
equation  (2.45)  can  be  written  as  a quadratic  form 


: (B-ift)  which  are 


expressed 


A can  be  expressed  as 


expressed  as  UTAU, 
of  three  matrices. 


no  reformulation  of  matrices  A2,A2,A2  is  needed  when  tl 
aspect  ratio  r is  changed.  This  permits  the  method  i 
consider  arbitrary  changes  in  aspect  ratios  without  ai 
additional  reformulation.  As  a result,  the  aspect  ratio  ci 

The  steps  involved  in  the  algorithm  are  given  i 
1.  Read  the  values  of  V number,  aspect  ratio  r,  and  ' 


t grid  points. 


Calculate 


and  U,  the  eigenvector  o£  B-1  A. 

6.  Calculate  b. 

7.  Print  the  values  o£  E(x',y')  at  the  grid  points 
normalized  propagation  constant). 


2.4.4  verification  and  Case  Studies 

The  algorithm  is  coded  in  FORTRAN.  The  main  program 
reads  the  index,  sets  up  the  A and  B matrices  and  calls 


separate  subroutines  for  matrix  inversion  and  eigenvalue 
estimation.  Double  precision  arithmetic  is  used  to  keep  the 
numerical  errors  in  the  matrix  inversion  as  small  as 

To  verify  the  accuracy  of  our  Variational-Finite 
difference  method,  a specific  asymmetrical  step  index 


channel  waveguide  <nchannel  = 
considered.  For  this  waveguide. 

respectively.  It  is  seen  that  results  are  in 

Finite  difference  method.  The  results  ol 
Finite  difference  method,  however,  are 
arbitrary  aspect  ratio  and  index  profile. 


ire  in  excellent 
the  Variational- 
the  Variational- 


The  cases  of  practical  importance  are  diffused  channel 
waveguides  with  index  profiles  functions  mentioned  below. 


a profile  in  equation  (2.54)  is  called  the  Gaussian- 

process  in  which  the  dopant  source  is  depleted  entirely 
during  diffusion.  This  condition  is  valid  for  indiffusion 

approximated  as  an  erfc-Gaussian  function  (equation 
r the  fundamental  model  are 


computed.  Figures  2.21  and  2.22  show  the  b-v 
characteristics  for  aspect  ratios  1 and  2 respectively.  In 
both  the  figures  it  is  seen  that  for  the  same  V number,  the 
step  index  waveguide  has  the  highest  b followed  by  the 
Gaussian-Gaussian  and  erfc-Gaussian  waveguides.  A similar 
trend  is  observed  in  the  planar  graded  index  waveguides  as 


illustrated  in  figure 
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Figure  2.11  Field  profiles  for  itF0,l,2  for  erfc  profile 


3Ali03dd3 


I 

I 

1 

i 

I 

1 


i 


1 


CHAPTER  THREE 
DIPFUSION  SIMULATION 


i -exchange 


The  study  of  ion-exchange/diffusion  in  glass  is  imp 
tant  because  this  phenomenon  controls  the  formation 
waveguides  on  the  surface  of  glass.  The  diffusion 
essentially  determines  the  device  parameters  sue 

controls  the  magnitude  of  the  index  change  at  the 
When  an  ordinary  soda-lime-silicate  glass  slide  is 
a molten  salt  bath  containing  monovalent  ions  such  as  Ag+, 
K+,  Li+  or  Tl+,  these  ions  exchange  with  the  sodium  ions 

the  melt  and  then  migrate  inside  the  glass  by  a process  of 
diffusion.  The  general  problem  of  diffusion  of  ions  in 


are  derived  as  follows. 

Consider  two  monovalent  ion  species  A and  8 which  are 
undergoing  inter-diffusion  in  glass.  Let  and  be  their 
be  their  self-diffusion  coefficients. 


equation  [He58] 


R " Gas  constant,  8.317  J/K-mole 


F = Faraday  Constant,  96500  J/V-mole 
T ~ Absolute  temperature, 


The  similar  equation  can  be  written  Eor 


Electroneutrality  requires  that  the  total  concentration 


and  (3.2),  an  expression  for  can  be  derived  as 


r-diffusion  coefficient 


diffusion  coefficient  DAB  is  concentration  dependent. 
During  the  fabrication  of  glass  waveguides  using  silver- 

bath  (and  hence,  in  the  waveguide)  is  much  smaller  than  that 
of  the  sodium  ions.  This  fact  considerably  simplifies  the 
diffusion  equation  as  follows. 


concentration  of  Ag' 
concentration  c 


As  mentioned  above,  CA  <<  CB,  with  the  result,  equation 
(3.6)  simplifies  as 


influence  c 


silver  self-diffusion 
inter-diffusion  process 


With  this  simplification,  several 
practical  fabrication  are  considered 
following  sections. 


as  a stronger 
t assumed  to  be 
ses  arising  in 


^-Dimensional  Diffusion  without 


i glas 


: diffusion,  the 


diffusion  can  thus  be  considered 


ion-exchange  takes 
ions  penetrating  to 


i-dimensional . The 


those  present  in  the  salt  bath,  so  that  the  concentration  c0 
of  silver  ions  in  the  melt  can  be  assumed  to  be  constant 


The  one-dimensional  diffusion  equation  with  constant 
diffusion  coefficient  is  (CrS61 


With  the  initial  condition,  c(x,0)  « 0 for  x > 0 and  the 
boundary  condition  c(0,t)  * cQ. 

closed  form  solution  of  this  equation  for  the  above 
mentioned  boundary  and  initial  conditions. 

Consider  the  function  c(x,t)  and  define  its  Laplace 


(3.12) 


t;  --'gt-  -5? -9 


r if « . i; . . r . 


CIO)  - /“  c(0,t)  e~3tdt  - /’  c0  e~3tdt  - |2 


a„  .£.-*sn 


ii  = mobility  of  silver  ions  in  glass 

and  mobility  vary  with  temperature  and  the  temperature 
dependence  is  given  by  an  Arhenius  type  relation  (Do641 


To  tailor  the  index  profile,  usually  the  electric  field 

to  resort  to  numerical  techniques.  A number  of  methods  have 

in  euqation  (3.22)  are  replaced  by  finite-difference 
approximations.  A rectangular  region  is  defined  in  the  t-x 
domain  with  t as  the  ordinate  and  x as  the  abscissa  as  shown 
in  figure  3.1.  Line  PS  represents  the  condition  x = 0,  so 
the  boundary  condition  c(0,t)  is  used  to  assign  the  values 


represents  the  condition  1 


(i,jj  are  then  approximated  as 


. c[  i+l.i)  - 2c(i,i)  ■■ 


Then  the  partial  differential  equation  (3.22)  is  transformed 
into  a set  of  coupled  difference  equations  by  usinq  the 
approximations  (3.25)  through  (3.27). 


c(i,j)  - e(i,J-l)  . D[c(i*lf j)  - 2cli.j)  * cd-lf.j), 


- c<i,j-l)  = <. 


* 6(j)  c(i-l,j) 


of  graphics  si 
and  surfaces. 


i program,  the  electric  field  is 

d diffusion  equation  (3.12)  for 
) form  analytical  solution  is  available.  The 


figure  3.2.  This  numerical  solution  shows  an  excellent 
agreement  with  the  analytical  solution  (complementary  error 
function)  establishing  the  validity  of  the  algorithm. 

The  program  is  then  used  to  calculate  the  concentration 
profile  resulting  from  different  values  of  E in  the  range  of 
0-100  volts/mm  and  t in  the  range  of  10-50  minutes.  These 


used  in  our  experiments.  Figures  3.3  and  3.4  show  the 


results  corresponding  to  various  process  conditions.  It  can 


time  or  a large  value  of  field  for  small  time  to  get 
approximately  the  same  concentration  profile.  To  be 
specific,  in  figures  3.3  and  3.4,  waveguides  with  E = 60 
v/mm  and  t = 30  minutes  have  approximately  the  same  profile 
as  those  with  E ■ 40  v/mm  and  t > SO  minutes. 


Figure  3.5  illustrates  how  the 
tailor  the  index  profile.  For  small 

function.  As  time  proceeds,  the 


values  of  time,  i.e.. 


profile  changes 


Gaussian-like  shape  and  subsequently  to  almost  the  step 
index  profile.  Thus,  the  simulation  has  shown  that  a 
judicious  combination  of  E and  t can  be  used  to  get  any 

profile,  however,  a large  value  of  field  should  be  used  for 
a large  diffusion  time. 


3.3  Two-Dimensional  Diffusion  With  And  without  Drift  Field 
In  the  previous  section,  one-dimensional  diffusion 

been  considered.  Channel  waveguides  in  glass  require 
diffusion  through  a window  which  requires  analysis  of  two- 
dimensional  diffusion.  The  diffusion  equation  for  the  two- 


and  the  region  near  the  top  surface  where  the  diffusion 
in  that  region  can  be  assumed  to  be  perpendicular  to  the 


simplif icatic 


The  window  size  has  a rather  small  ef£ect  on  the 
elective  depth  of  the  waveguide  for  both  ordinary 

sizes  w = 5 and  w = 10  pm,  the  effective  depth  is 
seen  in  figures  3.8  and  3.9  with  an  applied  field 


to  ■ 10  pm  for  the  same  t = 50  minutes. 

for  waveguides  fabricated  without  an  applied  field, 

them  quite  shallow.  Their  aspect  ratio,  which  is 

effective  depth,  is  also  • 2 or  more.  For  example 
in  figure  3.6,  effective  width  ■ 18  urn  and 

figure  3.7,  effective  width  • 10  urn,  effective 
depth  - 5 pm,  so  aspect  ratio  - 2. 
with  field  assisted  diffusion,  considerable  pene- 
tration of  silver  ions  without  sideways  diffusion 

In  figure  3.9,  w = 5 pm,  E = 30  v/mm,  t = 50  rain., 
effective  width  * 10  pm,  effective  depth  - 8 pm. 

an  optical  fiber  with  core  diameter  of  10  pm. 


results  indicate 


that  this  program  can  be 
effectively  used  to  design  channel  waveguides  exhibiting 
good  coupling  efficiencies  with  optical  fibers  of  specified 
core  diameters.  In  other  words,  the  values  of  mask  width  w. 


e resultant  waveguide 


illustrate  t 
mentioned  cas 

terminal.  The  visual  information  conveyed 
used  to  determine  the  effect  of  changing  tl 
and  t on  the  shape 


diffusion  t can  be  chosen  so  that 

of  the  fiber.  Figures  3.10-3.13 

running  the  program  for  different 
ie  figures  are  displayed  on  the  CRT 


facilitating  faster  design  procedures. 


Eectangular  grid  for  one  dimensional  dii 


mri  oi 


C(*.> 


Concentration  surface 


Concentration 


Concentration 


4 . 1 Fabrication  Process 

The  process  used  Cor  fabricating  ion-exchanged  glass 
waveguides  as  reported  in  literature  [Iz72,  St77,  Ch80]  is 
shown  in  figure  4.1.  It  has  a salt  bath  containing  a 

(NaN03)  and  10%  silver  nitrate  ( AgN03 ) heated  to  a 
temperature  of  320-360DC.  A soda-lime-silicate  glass 
substrate,  about  1 mm  thick,  is  suspended  in  the  bath  for  a 
period  of  40-60  minutes.  The  silver  ions  (Ag+)  in  the  bath 
replace  the  sodium  ions  (Na*)  in  the  glass  slide,  thus 

with  Ag+  exhibits  higher  refractive  index,  enabling  light 
guidance.  A variation  of  this  process  is  shown  in 
figure  4.2  in  which  a metallic  electrode  is  deposited  on  the 
back  surface  of  the  glass  slide  to  facilitate  the 
application  of  an  electric  field  which  enhances  the 
diffusion  of  silver  ions.  Standard  photolithographic 
masking  and  delineation  are  used  to  form  a window  or  strip 
pattern  using  aluminum  mask  as  shown  in  figure  4.2.  Ion- 
exchange  through  the  strip  region  is  used  to  fabricate 


several  sporadic  attempts  spread  over  about 
'82),  no  one  has  fabricated  low-loss  single 
dith  assured  reproducibility.  As  a result, 

a not  been  encouraging  from  a practical 
results  of  diffusion  simulation  and  propa- 
; of  this  project  indicate 


gation  analysis  performed  as  a pari 
that  the  silver  ion  concentration 
tant  parameter  1 

order  to  make  reproducible  waveguides, 
true  for  multimode  waveguides  and  bate) 

tions  concerning  the  existing  process  are  in  order: 


and  controlled  in 

processing  of  single 
i following  observa- 


2.  The  existing  process  takes  high  weight  ratio  of  AgNO^  to 
NaNO^  resulting  in  multimode  (7-16  modes)  waveguides. 

3.  With  such  a high  concentration,  reduction  of  silver  ions 

excessive  absorption  and  scattering  loss. 

4.  The  time  of  diffusion  necessary  to  produce  a single  mode 
waveguide  with  such  heavy  concentration  of  silver  ions 
is  only  about  2-3  minutes.  This  makes  the  timing  very 


making 


reproducibility 


5.  A single  mode  waveguide  fabricated  under  such  conditions 
is  usually  shallow.  Therefore,  such  a guide  couples 
poorly  to  an  optical  fiber  which  can  be  thought  of  as  a 
buried  structure  of  reasonably  large  cross-sectional 
dimension  (8-10  um) . 

necessary  to  start  out  with  very  small  concentration  of 
silver  ions  and  then  control  it  precisely  at  all  times 
during  diffusion.  The  reduction  of  Ag*  concentration  would 


60  minutes  even  for  a single  mode 
A novel  process  capable  of  releasing 


waveguides.  thus 
Paraday's  law  of 


electrolysis: 

When  96500  coulombs  of  charge  are  passed  through  a 

released  into  the  solution  as  ions. 

Thus,  generation  for  silver  ions  can  be  accomplished  by 
using  a metallic  silver  electrode  rather  than  a salt  mixture 

kept  minimal  as  ultra  pure  (99.999%)  silver  rods  are  commer- 
cially available.  The  process  of  releasing  silver  ions  can 
now  be  electronically  controlled  through  the  realization  of 


variable  current  source  using  off-the-shelf 
operational  amplifiers.  In  addition,  the  Ag*  concentration 
in  the  bath  can  be  changed  with  time  so  that  the  diffusion 

thus  providing  a larger  control  of  the  refractive  index 
profile.  The  basic  schematic  diagram  of  the  electrolytic 
process  is  shown  in  figure  4.3.  The  set-up  consists  of  a 
salt  bath  heated  by  a band  heater  coiled  around  it.  The 
temperature  of  the  bath  is  sensed  by  a thermocouple  and 
controlled  with  a closed-loop  controller  to  an  accuracy  of 
±0.5°C.  The  current  source  is  essentially  a voltage  to 
current  converter,  where  a Zener  stabilisied  input  voltage 
produces  precise  output  current  with  an  accuracy  of 
±0.1  mA.  The  redox  (reduction-oxidation)  reactions  at  the 
electrodes  are 


Anode  (silver  rod) 


Cathode  (platinum  wire) 


9 (metal) 
a(metal) 


since  sodium  is 
to  be  deposited  as 


fritted  glass  enclosure  is  used  around  t 
shown  in  figure  4.4.  The  bore  size  if 
is  " 2-4  im.  It  essentially  prevents 

sodium  ions  are  made  to  react  at  the 
sodium  which  collects  inside 


:he  platinum  wire  as 
n the  fritted  glass 

! cathode  producing 


enclosure  [Na84] . The  rate 

which  they  are  generated  by  electrolysis,  s 

This  novel  technique  has  an  inherent  ability  to  generate 

capability  is  extensively  used  in  this  project  work  to 
facilitate  fabrication  of  reproducible  single  mode 
waveguides. 


Concentration 


The  electrolytic  process  which  can  generate  precise 

on-line.  The  conventional  methods  to  measure  concentration, 
e.g.,  mass  spectrometry,  atomic  absorption, 

suitable  as  they  are  off-line,  thus 


relatively  low  temperature  (0-80"CJ  and  therefore  are  useful 
gh  temperatures  in 


situ  has  been  developed  [La84c].  The  operation  of  the 
transducer  is  described  in  the  remainder  of  this  section. 

As  shown  in  figure  4.5.  the  glass  bulb  transducer  is 
basically  a galvanic  cell  in  which  the  redox  reaction  is 


3 particular  cell  the  s 
anode  and  cathode  compartments. 


The  generated  cell  potential.  called 
potential,  is  given  by  [B159] 


reference  it 


P substituted,  equation  (4.1)  c 


Concentration  cells  u 
: glass  bulb  or  glass  t 
ive  been  reported.  It 


g fritted  glass  [Du61]  in  place 
s with  asbestos  fibers  (F1561 
our  experience,  however,  that 


fritted  glass,  complete  mixing  between  the  solutions 
cell  potential  monotonically 

the  glass  bulb  itself  can  be  used  as  a 
l [In621  has  reported  the  use 
a containing  reference  electrode  to  form 
The  underlying  idea  is 
: elevated  temperature  ( 330°C-370°C)  can 
ion-bridge  between  the  solutions  since 
s rather  high  mobility  at  these  temper- 
•e  no  pores  present  within 
e physically  isolated. 

i Pyrex  glass  tes 
approximately 


form  a glass  bulb  with  a tubular  stem.  The  glass  bulb  is 
nearly  spherical  in  shape  with  about  12  mm  diameter  and 
about  0.1  mm  wall  thickness  as  shown  in  figure  4.6.  A 
reference  melt  is  produced  by  releasing  2.0  gms  of  silver 
from  a high  purity  (99.999%)  rod  into  a melt  of  400  gms 

is  approximately  equal  to  3.94  x 10"3  mole  fraction  (MF). 

4 gms  of  the  reference  melt.  High  purity  silver  wires  are 

the  anode  immersed  in  the  bath.  The  cell  is  then  used  to 


impedance  digital  voltmeter. 

By  electrolytically  releasing  known  quantities  of 
Ag*  ions  in  400  gms  of  NaN03,  different  known  concentrations 

potential  each  time,  a plot  of  cell  emf  vs.  log(c3/C2),  as 
shown  in  figure  4.7  is  obtained.  The  solid  line  represents 
the  theoretical  cell  potential  as  given  by  equation  (4.2). 
As  expected,  both  curves  exhibit  linear  behavior.  Both  the 
experimental  and  theoretical  lines  have  nearly  the  same 
slope,  with  the  slope  of  the  experimental  line  slightly 

direction  amounts  to  almost  20°C.  The  bath  temperature  is 


between  the  contents  o£  the  glass  bulb  and  that  of  the  bath, 
causing  the  difference  in  the  slope  of  the  experimental  and 
the  theoretical  line.  Moreover,  when  two  concentrations  are 
equal,  an  asymmetry  potential  of  about  14  mV  is  measured 

[En621  attributes  the  asymmetry  potential  to  a difference  of 
sodium  ion  concentration  in  the  glass  between  the  inside  and 
outside  of  the  bulb  which  occurs  during  blowing,  we  believe 


3 possible  to  implement 


amplifier  is  i 
configuration,  th 

already  in  place, 


.e  concentration  transducer, 

log  term  in  it,  so  that  an  anti-log 
jet  a signal  which  varies  linearly 
In  the  feedback 

w the  preset  value,  additional  silver  ions 
by  periodically  pumping  currrent  in  the 

controls,  some 


tation  of  the  entire  process  control  will  mak 
to  fabricate  glass  waveguides  reproducibly  ! 


: possible 


process  variables 


controlled, 


.3  Planar  Waveguide  FabcicaCion  Procedure 
silver  ion  releasing 


de-ionized  (DI)  water.  Remove  the  surface  layer 
inside  wall  of  the  beaker  by  boiling  dilute 
etch  (16-1-1-2  etch):  16  parts  phospheric  acii 

acetic  acid,  1 part  nitric  acid  and  2 parts 
for  10  minutes.  Rinse,  with  DI  water.  Dry  in 
at  150°C  for  5-6  hours  (preferably  overnight). 

( A.C.S-grade,  Fisher)  using  handgloves.  Take 

beaker  to  330°C.  Handling  with  gloves,  weigh  t 
rod  and  note  its  initial  weight.  Add  pure  NaN 
fritted  glass  tube  and  insert  the  platinum  wi 


the  tube.  Suspend  the  silver  rod  and  the  fritted  glass 
tube  in  the  molten  NaNO^  and  wait  till  thermal 
equilibrium  is  reached.  By  pumping  current  through  the 

for  60  minutes  releases  0.4  gms  of  silver  in  the  bath. 

release  any  arbitrary  amount  of  silver  ions  in  the 
melt.  During  the  ion  release,  no  precipitation  should 
be  observed  and  the  molten  salt  should  be  clear.  If 
not,  some  silver  is  being  deposited  back  to  metallic 


glass  disk.  After  the  time  has  elapsed,  reduce  the 
current  to  zero,  remove  the  silver  rod  from  the  melt, 
rinse  it  with  hot  DZ  water  and  weigh  it  again.  Compare 
the  weight  loss  with  the  theoretically  calculated  weight 


(B)  Diffusion 

Clean  the  glass  slides:  Boil  in  trichloroethylene 

(TCE),  acetone  and  methanol  for  10  minutes  each.  Rinse 
in  warm  DX  water  and  spin  dry.  Hardbake  at  150aC  for 

60  minutes. 


of  diffusion. 

Remove  the  slides  from  the  bath  after 


the  required  time 
diffusion  and  let 


them  cool  for  10  minutes. 

Remove  the  solidified  salt 
waveguides)  by  shaking  them  in  w 
waveguides  with  soft  tissue  paper 
Clean  the  waveguides  surface 

of  the  mode  indices. 


coupling  light 


Existing  process  for  glass  waveguide 
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Electrolytic 


releasing 


Ag 
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4.4  Introduction  of  fritted  glass  around 
cathode  in  the  electrolytic  process 


potential 
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CHAPTER  FIVE 

PLANAR  WAVEGUIDE  CHARACTERIZATION 


evaluation  c 
refractive  i 


most  useful  techniques  employed  in  the 
s propagation  constant  of  planar  surface 
: prism  coupling  technique.  As  shown  in 


wavelength  of  the  incident  light.  > 
energy  tunnels  through  the  gap  region 
guiding  mode.  Thus*  each  guiding  mode  i 

angles  enables  one  to  calculate 


1 compared  to  the 

d couples  to  the 

measurement  of  these 
effective  mode  index 


incidence 


refraction 
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The  travelling  wave  in  the  prism  has  a wave  vector  with 
magnitude  k0n-  which  can  be  decomposed  into  a horizontal 
component  kQnD  cos  0 and  a vertical  component  kQn_  sin  0. 
The  tunneling  interaction  takes  place  when  the  horizontal 
component  o£  the  wave  vector  (along  the  z direction)  is 
equal  to  the  propagation  constant  of  a given  mode  among  the 
many  possible  modes  of  the  waveguide.  Under  this  condition 
that  particular  mode  gets  excited  (Ti71). 


of  incidence  and  the  prism 


a precision  rotating  table  as  shown  in 

figure  5.1  is  used. 

The  angle  measurement  is  carried  out  w 

ith  a resolution  of 

6 arc  seconds.  The  prism  refractive  ind 

ex  is  measured  to  an 

accuracy  of  4 decimal  places  using  an  . 
For  glass  waveguides  (substrate  refract 
a Schott  glass  (SF18)  prism  with  n. 

Abbe'  refractometer. 

0.6328  urn)  is  used.  An  output  prism  ca 

in  be  used  to  couple 

sensitivity  of  the  method,  the  beam  is  mechanically  chopped 


ation  is  used  in  conjunction  with  a 

1 precision  lock-in 

amplifier  (Princeton  Applied  Research  He 

Jdel  124A)  to  detect 

the  guided  light  amplitude.  The  point 

of  maximum  coupling 

is  determined  as  the  input  incident  angle  is  slowly 


increased  by  rotating 


avoiding  any  backlash 


may  be  present  in  the  rotating  table.  This  method  greatly 


propagation  loss  in  optical  waveguides:  absorption  losses 

and  scattering  losses. 


tal  and  hydroxyl 
y controlling  the 
i-clean  processing 

1 impurities  less 


ions  (0-H).  These  losses  are  kept  small 
composition  oC  the  substrate  glass.  Ult 
environment  is  used  while  fabricating  t 

(e.g.,  Schott  S-8011)  w 
than  a few  parts  per  million. 

The  scattering  losses  again  are  of  two  types:  volume 

etc.,  within  the  volume  of  the  waveguide.  This  type  of  loss 
depends  very  strongly  on  the  relative  size  of  the  imperfec- 
tions in  comparison  with  the  wavelength  of  light  in  the 

as  few  in  number.  As  a result,  the  volume  scattering  loss 
is  negligible  compared  to  the  surface  scattering  loss. 


Surface  scattering 
for  the  fundamental  me 
waveguide  because  most 


confined  to  the 
happens  to  be  ne 

fabrication  becau 
waveguide  should 


be  significant  particularly 

region  of  higher  refractive  index  which 

important  in  the  waveguide  design  and 
*e  from  an  application  point  of  view,  the 


waveguides, 
in  which  a m 


propagation  loss 
in  pcinciple,  u 


in  planar  surface 


distance  hai 
absorbed  by 
magnitude  of  the 

necessary  to  keep 

rather  restrictive 
while  moving  t 


by  moving  the  output  prism  along 
waveguide.  The  plot  of  output  power 

e impurities  or 

ep  the  same  coupling 
output  prism 
input  prism  t 


scattering. 


condition  seems 
> get  disturbed 


n addition,  extreme  c 


e output  prism, 
ided  to  avoid  spoiling  o 

An  improved  3-prism  method  (W08O]  in  which  the  results 
? independent  of  the  prism  coupling  coefficients  is  used 


Let  Pa  and  be  the  output  powers  of  prisms  A an* 
respectively.  ya  and  Ya  are  the  coupling  coefficients 
the  output  prisms. 


e attenuation  in  th 


P(S2)  = IP(Zj) 


output 


rBIPUi) 


4pb  ’ pb  - pbo  ■ Va  p'*i> 


Optical  interferometric  technique:  In  this  method,  the 

index  profile  is  determined  by  measuring  the  optical 
phase  shift  induced  by  the  refractive  index  variation. 
The  method  uses  a Mach-2ehnder  or  a similar  type  of 
interferometer  to  examine  the  polished  slice  of  the 
waveguide.  This  technique  is  applicable  essentially  to 
relatively  deeper  waveguides  (few  tens  of  microns), 

effective  depths  in  the  range  of  6-10  urn,  cross-section 
preparation  is  somewhat  difficult  because  of  the 

polishing  thin  slices  of  the  waveguides. 

Surface  analytical  techniques:  Surface  analysis 

involves  the  determination  of  silver  ion  concentration 
profile  by  using  one  of  the  destructive  surface  analysis 
methods,  e.g..  X-ray  microprobe,  Auger  spectroscopy, 
Secondary  ion  mass  spectroscopy,  etc.  The  implicit 

concentration  profile  matches  the  refractive  index 

exchanged  waveguides  because  the  presence  of  silver  ions 
is  solely  responsible  for  an  increase  in  the  refractive 


project  deals 


depths  are  typically  in  the  range  of  6-8  urn,  the  optical 
interferometric  techniques  were  not  used  for  index  profile 
determination.  A Commercial  Mach  Zehnder  interferometer  is 
also  not  available.  Several  other  groups  working  on  glass 
waveguides  have  been  using  Electron  Probe  Micro  Analysis 

(St77 , LiB2] . The  equipment  needed  for  this  technique  is 
not  available  at  the  university  of  Florida.  A new  and 
simple  method  using  scanning  electron  microscope  (SEM)  is 
developed  for  the  index  profile  measurements. 

In  an  SEM,  the  specimen  is  bombarded  with  an 
accelerated  beam  of  electrons  with  the  accelerating  voltage 
in  the  range  of  20-50  KV.  The  specimen  generates  two  types 


electrons  [Go81] . 

The  secondary  electrons  are  low  energy  electrons  { less 
than  50  eV),  generated  primarily  by  inelastic  scattering  of 
the  incident  electron  beam.  Since  the  secondary  electrons 
are  generated  with  low  kinetic  energy,  they  usually  come 
from  the  surface  and  carry  information  about  the  surface 
topography.  Their  yield  is  relatively  insensitive  to  the 

with  the  atomic  number  of  specimen. 


The  process  of  back  scattering  usually  takes  place  as  a 
result  of  a sequence  of  elastic  scattering  events  in  which 
the  net  charge  in  the  direction  is  sufficient  to  carry  the 
incident  electrons  out  of  the  specimen.  The  back  scattered 
electrons  ace  high  energy  electrons  {greater  than  100  eV) 

electron  coefficient,  defined  as 


nBs  = number  of  backscattered  electrons 

rig  b number  of  beam  electrons 

iBs  B back  scattered  electron  current 


If  the  backscattering  coefficient  n is  measured  as  a 


(5.16) 


when  the  target  is  a homogenous  mixture  of  several 
elements,  a simple  rule  of  mixture  based  on  the  weight 
fraction  of  the  components  is  found  to  apply 


element  and  is  the  weight  fraction  of  that  element. 
Experimental  measurements  indicate  that  the  backscattering 
coefficient  does  not  depend  upon  the  beam  energy  [Go8 1 ] . 

In  the  Ag+  - Na+  exchanged  waveguides,  the  substrate 
material  mainly  contains  ions,  such  as  Na+,  Si4+,  0 , Ca++, 

(Ag+).  The  atomic  number  of  Ag  (47)  being  higher  than  that 
of  Na  (12)  and  Si  (14),  the  back-scattering  coefficient  of 


silver  ions  in  glass  (which  is  true  in  the  case  of  glass 
waveguides),  still  contributes  significantly  to  the  yield  of 
the  backscattered  electrons.  Thus,  if  the  sample  is 
polished  across  the  edge  of  the  waveguide  and  scanned  in  the 
SEM,  the  backscattered  electron  intensity  profile  gives  a 
good  approximation  of  the  silver  ion  concentration  profile 
and  thus,  of  the  index  profile.  The  SEM  is  calibrated  with 


and  the  diffusion  depth 


The  procedure  can  be  outlined  as  follows.  The 
waveguide  is  polished  in  a plane  perpendicular  to  the 
guiding  direction.  A thin  layer  of  epoxy  is  applied  to  the 
surface  of  the  waveguide  and  a glass  slide  made  from  the 

stuck  to  the  surface  before  polishing.  This  prevents 
rounding  or  chipping  at  the  edges  of  the  waveguide  while 
polishing.  The  polishing  is  done  in  two  steps:  Rough 

diamond  paste.  Glass,  being  a dielectric  material,  needs 
some  conductive  coating  to  drain  the  charge  built  up  from 
the  incident  electron  beam.  A carbon  coating  (50  A thick), 
has  produced  good  results.  Since  the  atomic  number  of 
carbon  (6)  is  much  lower  than  that  of  silver  (47),  the 


carbon  coating  does  not  interfere  with  the  measurements. 

electron  and  the  backscattered  signals  are  displayed 
simultaneously  on  the  output  screen  split  in  two  halves.  In 

for  composition  and  Topo  for  topography.  Figure  5.6 


a "Compo  + Topo"  photograph.  The  left  half 
» the  secondary  electron  signal  or  "Topo" . 
any  pattern  indicates  that  there  is  no 

right  half  of  the  picture  is  the  display  of 


presence 
vertical 
high  molecular 


silver  ions  diffused  in  the  glass, 
e represents  the  epoxy  layer.  The  e 


little  to  the  yield  of  backscattering. 
of  the  backscattering  signal  gives  th 
shown  in  figure  S.7.  The  big  dip  in  thi 


3 signal  corresponds 

the  surface  of  the  waveguide.  Such  a reference  point  is 

The  profile  in  figure  S.7  is  of  a waveguide  fabricated  by 
silver  ion  diffusion,  without  an  applied  electric  field.  It 
looks  essentially  like  a complementary  error  function, 

diffusion  simulation,  as  reported  in  Chapter  Two.  The 
profile  is  calibrated  in  the  depth  direction,  i.e.,  a 
L line  represents  a distance  of  10  urn  as  shown  in 


equation  3 
deeper  (~  : 


simulation, 


graph.  Figure  5.8 
essentially  like  a step 


i diffusion.  The 
gain  confirms  the 


with  field  assisted  diffusion  are  considerably 
22  urn),  than  those  made  with  ordinary  diffusion 
This  again  confirms  the  results  of  diffusion 
, as  can  be  seen  by  comparing  figures  3.2  and  3.3. 
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T'opo"  photograph  bj 


scanning 


fabricated  by 


diffusion  without  applied  field 


CHAPTER  SIX 

PLANAR  WAVEGUIDE  DESIGN  PROCEDURE 

fi.l  Relations  Between  the  Process  and  the  Device  Parameters 
In  this  section,  experimental  results  of  the  diffusion 
studies  and  modal  behavior  of  planar  glass  waveguides  is 
described.  The  planar  glass  waveguides  fabricated  by 
diffusion  without  an  applied  electric  field  have  a 

theoretical  treatment  presented  in  Chapter  Three 
(figure  3.2)  has  been  confirmed  experimentally  using  SEN 

characteristics  for  this  highly  asummetric  profile  have  been 
determined  by  solving  the  normalized  mode  dispersion 
equation  in  Chapter  Two  (figure  2.9).  A combination  of 
these  theoretical  results  and  experimental  investigation  is 
used  to  find  the  links  between  the  process  and  the  device 
parameters.  The  experimental  procedure  is  outlined  below. 

Several  sets  of  single  mode  and  multimode  waveguides 
are  fabriated  using  Labmate  microscope  slides  for  different 
values  of  the  process  parameters,  CQ  (Ag+  concentration  in 
bath)  and  t (diffusion  time).  For  a particular  batch,  a 
known  value  of  Cc  is  generated  using  the  electrolytic 


fabricated 


coupling  technique  described  in  Chapter  Five  is  used  to 
excite  the  modes  and  to  estimate  the  mode  indices  (equation 
(5.8)).  A computer  program  is  written  to  estimate  the 
diffusion  constant  D and  the  maximum  index  change  An.  The 
procedure  is  described  as  follows. 

As  before,  the  V parameter  is  defined  as 


= waveguide  depth 

= maximum  change  in  the  refractive  i 


■ function  t 


complementary 


wavguide 


diffusion  coefficient 


t « time  of  diffusion. 


measured  up  t 
refractometer. 


urn.  The  bulk  refractive  inc 
* Labmate  glass.  nb  = 1.5125. 


Now  the  normalized  mode  index  b for  an  asymmetrical  graded 
index  waveguide  is  defined  as 


surface  ii 


approximation  {which 


From  equations  (6.4)  and  {6.6}  it  is  seen  that  Cor  a 
given  glass  substrate  and  a given  wavelength  oC  operation,  b 
and  V are  functions  oC  D , An  and  t.  An  depends  upon  C_,  the 
concentration  oC  Ag+  in  bath.  Thus,  a set  oC  waveguides 


fabricated  with  the  same  CQ 
same  An  but  different  depths, 
The  computer  program  reads  as 
mode  indices  and  diffusion  tit 


fitting  experimental  data 
characteristics  (figure  2.9).  The  val 
initially  assumed  and  for  this  chc 
calculated  from  equations  (6.4)  and  ( 
defined  as  the  largest  difference  betwet 

then  varied  over  a prespecified  r« 
calculations  are  repeated.  The  errc 
points  in  the  An-D  plane  are  plotted 
6.1.  The  contours  show  that  the  error  1 


or  algorithm  is  used  for 
theoretically  derived  b-V 


s calculated 
f calculated 


and  D corresponding 


The  waveguides  are  fabricated  with  several  values  of 
C-,  and  an  is  estimated  for  each  CQ.  The  relationship 
between  an  and  CQ  is  plotted  as  shown  in  figure  6.2.  The 

(1-6),  the  relationship  between  an  and  CQ  can  be 
approximated  as 


i is  an  important  linl 
ie  process  parameter  C 
r Labmate  slides  that 


ween  the  device  parameter 
:his  empirical  relation  is 
used  in  our  experiments. 
3 weight  fraction  of  Na^O 


in  the  substrate  material  is  changed.  However,  for  the 
weight  fraction  in  the  neighborhood  of  10-15%,  an  order  of 

To  estimate  the  variation  of  the  diffusion  coefficient 

fabricated  at  330,  354  and  364°C.  The  diffusion  coefficient 

mentioned  above.  The  variation  of  D with  T is  given  by  an 
e relation  [Do641 , 


Energy  and 


Taking  logarithm 


diffusion  time  is  of  the  order  of 
diffusion  depth  is  of  the  order  of  a few  microns 
expressed  in  micrometervminute.  The  temperature  r 
restricted  as  NaN03  melts  at  307°C  and  decomposes  a 
[CR831.  Figure  6.4  shows  the  plot  of  D vs.  T o- 
relevant  temperature  range.  It  is  seen  that  the  cu 
be  approximated  by  two  linear  segments  with  a trs 


Beyond 


that  in  this  region,  the  temperature  fluctuations  during 
fabrication  will  cause  significant  changes  in  D resulting  in 
a large  variation  in  the  guide  depth  d.  This  indicates  that 
the  diffusion  temperature  should  be  chosen  between  307®C  and 
353®C.  During  the  experiments  it 
temperature  gradient  of  as  much  as  2( 


diffusion  temperature,  which  is  sufficiently  higher  than 
307®C  (the  melting  point  of  NaNOj). 

At  330*0,  the  diffusion  coefficient  has  a value. 


Substituting 


This  is  the  second  important  link  between  the  device 


6.2  Design  Procedure 

Using  the  relations  6.7  and  6.10,  a systematic  design 
procedure  is  formulated  as  described  below  [La84d] . 

The  normalized  dispersion  equation  for  an  asymmetrical, 
graded  index  planar  waveguide  is  given  by  equation  (2.11). 


function  pcofilej 


Thus,  Vcm,  the  value  of  V number  at  cut-off  of  the 
mode  is  given  by 


vcm  f /erfc(x') 


since  erfc(x')  is  a monotone  decreasing  function 
erfc( 2)  » 0.0005, 


From  equation  (6.3), 


from  equations  (6.13) 


(6.14), 


.-4S&- 


.7)  is 


experimentally  confirmed  by  fabricating  various  graded 
(erfc)  waveguides  by  electrolytic  release  technique 


prism  coupler. 


Ag+  MOLAR  CONCENTRATION 


Figure  6.5  Design 
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description 
study  of  Ag4 


project  is  described,  followed  by  a brief 
the  research  that  is  needed  to  complete  the 
a+  ion-exchanged  waveguides. 


7 . 1 Mode  Profile  Measurement 
Mode  profile  is  an  important  propagation  parameter  for 
both  the  surface  and  buried  waveguides.  When  two  waveguides 
are  connected  end  to  end  for  coupling,  the  efficiency  of 
power  transfer  depends  upon  the  overlap  of  the  fields  in 
these  waveguides.  Thus,  in  a variety  of  coupling  problems 
such  as  laser  to  fiber,  laser  to  waveguide,  fiber  to 


waveguide  and  waveguide  to  waveguide  coupling,  it  is 
mode  profiles.  Two  methods  are  typically  used  to  carry  out 
short  focal  length  placed  near  the  waveguide  to  magnify 


plane  and  an  oscilloscope  to  view  the  video  signal  [Xo82]  , 
a suitable  semiconductor 


employs 


manually  [Ke79] . 

The  latter  technique 
manual r hence  slow.  Wit 


small  aperture  t 


disadvantage  of  being 
be  viewed  continuously 


profiles  i 


to  smooth  out  the  electronic 
be.  This  technique  is  used 
planar  surface  waveguides 
set-up  is  shown  in  figure  7,1. 
propagating  three  modes  is  polished  at 
avoid  scattering  and  is  excited  using  a prism  coupler 
input  end.  The  modes  are  individually  excited  a 


observe  the  mode 
The  experimental 

;he  output  end  to 


waveform  calculation  capabilities  is  used  to  smooth  the  data 


profiles  for  the  first  three  modes.  Since  the  waveguide  is 
made  with  silver  ion  diffusion  without  drift  field,  it 
essentially  has  a complementary  error  function  type  index 

Chapter  Two.  As  seen  in  figure  2.11,  for  the  higher  order 
modes,  the  inner  lobes  of  the  field  profile  have  higher  peak 
amplitudes.  This  theoretical  result  is  confirmed 
experimentally  as  seen  in 


This  video  camera-oscilloscope  combination  can  be  used 
to  determine  the  mode  profiles  of  channel  waveguides  and  to 
measure  their  effective  modal  widths.  These  data  are 


from  buried  waveguides  [Ch80,  Li82) . In  this  process,  a 
surface  waveguide  is  formed  by  a field  assisted  diffusion  of 

field  assisted  diffusion  in  pure  NaNO^ . The  buried 
waveguide  is  formed  as  sodium  ions  replace  silver  ions  at 
the  surface.  The  following  four  process  parameters  control 


t^  : Time  of  diffusion  in  the  first  step  ( NaNO^  with  Ag+) 
: Field  applied  in  the  first  step 
t2  : Time  of  diffusion  in  the  second  step  (pure  NaNO^ ) 

&2  : Field  applied  in  the  second  step. 


In  the  existing  process,  field  in  the  range  of 
60-100  v/mm  is  used.  This  results  in  a step-like  waveguide 
(as  shown  in  figure  5.8),  which  when  buried,  produces  a 


waveguide 


found  experimentally  that  if  no  field  is  applied  in  the 

symmetric,  fiber-like  index  profile  as  shown  in  figure 
7.6.  This  is  the  first  time  a fiber-like  profile  has  been 
achieved  in  a planar  substrate.  Further  studies  are 
neaessary  to  understand  the  diffusion  process  used  for 


7.3  Channel  Waveguide  Fabrication 
For  the  fabrication  of  channel  waveguides,  the  same 
electrolytic  process  is  used  for  diffusion  in  which  the 
silver  ions  are  allowed  to  diffuse  only  through  the  windows 

pattern  fabrication.  Some  of  the  older  techniques  such  as 
wet  chemical  etching  which  have  poor  resolution  have  been 
replaced  by  a superior  lift-off  technique  which  is 

material,  but  other  metals  such  as  Titanium  could  also  be 
used.  Figure  7.8  shows  the  SEM  picture  of  a pattern  made  by 
lift-off.  The  high  quality  of  the  pattern  can  be  clearly 
seen  from  the  straightness  of  the  windows  and  the  sharpness 

due  to  lack  of  extremely  clean  photolithographic 
facilities.  It  is  also  necessary  to  determine  the  optimum 


combination  of  the  process  parameters  {photoresist  spinning 
speed,  photoresist  spinning  time,  exposure  time  and 
development  time)  to  improve  the  yield.  The  substrates  must 

this  cleanliness  should  be  maintained  throughout  the 
processing.  Any  dirt  or  micron  size  dust  particle  causes 

into  broken  photoresist  patterns. 


7.4  Future  Applications 

The  knowhow  generated  in  this  project  can  be  utilized 

insertion  loss  figure  of  1 db  for  1.3  cm  waveguide  [Ra82] 
can  be  improved  by  fabricating  high  quality  single  mode 
channel  waveguides  on  glass  substrate.  Integrated  optical 

couplers,  wavelength  filters,  etc.  can  be  designed  and 
fabricated  after  perfecting  the  techniques  for  channel 
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Figure  7.2  ttode  profile 
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PHOTOLITHOGRAPHIC  LIFT-OFF  TECHNIQUE 


Aluminum  deposition:  Deposit  about  500- >000  ar.gstorms  of  Al. 


Lift-off:  Use  warm  acetone  and  ultrasonic  vibrations  to  dissolve 
the  photoresist  patterns. 


Figure  7.7  Photolithographic  lift-off  technique 


Aluminum 


pattern  for  channel  waveguide  fabricatic 


CHAPTER  EIGHT 
CONCLUSIONS 


The  aim  of  this  project  was  to  achieve  a basic 
understanding  of  diffusion  and  optical  wave  propagation  in 
glass  and  to  use  the  results  of  this  research  to  fabricate 
single  mode  waveguides  with  assured  reproducibility.  A 

and  experimental  investigation  is  used  to  analyze  silver  ion 
diffusion  and  mode  propagation  in  optical  glass  waveguides. 

systematic  procedure  for  designing  single  mode  waveguides 
entirely  in  terms  of  the  process  parameters. 

: is  given  below: 


1.  Normalized  propagation  constants  and  mode  profiles  for 
arbitrary,  highly  asymmetrical  graded  index  planar 
waveguides  along  with  a closed  form  expression  for  the 
b-V  characteristics  for  well-guided  modes  in  buried 

2.  Propagation  analysis  of  channel  waveguides  with 
arbitrary  index  profiles  and  aspect  ratios  using  a new 


Variational-Finite  difference  method  resulting  in 
normalized  solutions  for  the  first  time. 

Diffusion  simulation  with  and  without  applied  electric 
field,  for  one  dimensional  and  two  dimensional  diffusion 
of  ionic  species,  e.g.,  Ag+  ions  {very  low 
concentration)  in  a binary  exchange  scheme  with  Na*  ions 

Design  of  a novel  electrolytic  process  for  releasing 
silver  ions  in  the  molten  NaNO^  bath  and  development  of 
a Galvanic  cell  transducer  for  in  situ  measurement  of 


Fabrication  of  planar  single  mode  waveguides  with  very 
low  propagation  loss — 0.2  db/cm  and  experimental 
confirmation  of  mode  profiles  in  graded,  index  planar 
waveguides. 

Estimation  of  the  diffusion  coefficient  and  empirical 
derivation  of  an — CQ  curve  for  Labmate  glass  (Height  % 


procedure  for  designing  single  mode  planar  waveguides. 


diffusion  process. 

Index  profile  measurement  by  a new  and  simple  technique- 


-scanning  electron  microscopy. 
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